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ABSTRACT 

We analyze the X-ray scattering halos around three Galactic Anomalous X-ray Pulsars (AXPs) 
in order to constrain the distance and the optical extinction of each source. We obtain surface 
brightness distributions from EPIC-pn data obtained with XMM-Newton, compare the profiles of 
different sources, and fit them with a model based on the standard theory of X-ray scattering by 
dust grains, both for a uniform distribution of dust along the line of sight, and for dust distributions 
constrained by previous measurements. Somewhat surprisingly, we find that for all three sources, 
the uniform distribution reproduces the observed surface brightness as well as or better than the 
distributions that are informed by previous constraints. Nevertheless, the inferred total dust columns 
are robust, and serve to confirm that previous measurements based on interstellar edges in high- 
resolution X-ray spectra and on modelling of broad-band X-ray spectra were reliable. Specifically, 
we find A v ~ 4, 6, and 8 mag for 4U 0142+61, IE 1048.1-5937, and 1RXS J170849.0-400910, 
respectively. For IE 1048.1—5937, this is well in excess of the extinction expected towards a H I 
bubble along the line of sight, thus casting further doubt on the suggested association with the 
source. 

Subject headings: stars: individual (4U 0142+61, IE 1048.1-5937, 1RXS J170849.0-400910) - stars: 
neutron - dust, extinction - scattering 



1. INTRODUCTION 

Among the diversity of objects comprising the neutron 
star family, Anomalous X-ray Pulsars (AXPs) are one of 
the two classes of magnetars, with their emission pow- 
ered by the decay of extremely strong magnetic fields, 
of 10 l4 -10 16 G, whose origin is still a m atter of debate 
(for a review. IWoods fe Thompson! [2006). As remnants 
of massive, short-lived stars, these sources lie preferen- 
tially near the Galactic Plane and thus are strongly ab- 
sorbed, hindering measurements of their intrinsic spec- 
tra and luminosities. In some cases, distances and ex- 
tinction columns could be estimate d from associated su- 
pernova remnants (IE 2259+58 6. iGregorv fe Fahlmanl 
19801: IKothes. Uvaniker fe Yarl |2002J; IE 1841-045, 
Vasisht fc GotthelJ I1997T and o ther structures (IE 
1048.1-5937. IGaensler etall[2005h. 

M ore recent ly, iDurant fe van Kerkwijk! (2006a, bj, here- 
after DvK06a,b) tried to determine the extinctions and 
distances to the sources more directly, by measuring the 
extinction from absorption edges in high-resolution X- 
ray spectra, and inferring distances by comparison with 
the run of extinction with distance as inferred from red 
clump stars near the line of sight. They found that all 
AXPs were located in spiral arms (not surprising given 
their youth) and that their luminosities were remark- 
ably uniform, near 10 ergs -1 . Comparing the inferred 
distances to earlier work, the results were consistent for 
IE 1841-045 and XTE J1810-197 (for the latter a dis- 
tance from H I absorption was available), but incon- 
sistent for IE 2259+ 586 and IE 1048.1-5937. For IE 
2259+586, DvK06a argued that the associated super- 
nova remnant CTB 109 could be at the larg er distance of 
~ 7.5 k pc instead of the ~ 3 kpc inferred by IKothes et al.l 
(|200H . while for IE 1048.1-45937, they concluded that, 

rivera,mhvk@astro. utoronto. ca 



at ~ 9 kpc, the source was not associated with th e H I 
bubble at ^2.7 kpc found by Gacnsle r et al.l ([2005D . even 
though it appeared fairly nicely centered on it. 

Of course, it is possible that the analysis of lDvK06al b 
was flawed. The underlying assumptions are that the 
extinction measures are reliable, that they reflect only 
interstellar extinction, and that the run of interstellar ex- 
tinction with distance can be reliably inferred from red 
clump stars. The latter assumption is perhaps the safest, 
since the red clump stars are fairly well-understood stan- 
dard candles. The extinction measures are relatively un- 
certain, however, especially for IE 1048.1—5937, which 
relied on fitting of the broad-band X-ray spectrum and 
is thus s ensitive to the assumed intrinsic model (see 
lDvK06bl ). Furthermore, if a large part of the extinc- 
tion were local to the source, and thus not present along 
the lines of sight to the red clump stars, the distances 
would be overestimated. 

Here, we seek to test the extinctions and distances 
using X-ray scattering halos, which appear as dif- 
fuse emission around an X-ray point source and arise 
from small-angle scattering of sof t X-rays by dust in 
the line of sight (|Overbeckl 119651 ). X-ray halo anal- 
yses have been carried o ut on a variety of galactic 
and extragalactic object s (|Mauche fe Gorensteinl 119861 : 
iPredehl fe Schmittlll995l ). and have proven to be partic- 
ularly useful in constraining the properties of dust grains 
in the ISM. For variable sources, where the time delays 
in the scattering halo could be mea sured, they have also 
been used to determine di stances (jPredehl et al.1 l2000l 
IXiang. Lee, fc Nowaldl2007l) . 

For our purposes, the benefit of the scattering halo is 
mostly that it provides a column density that is purely 
interstellar, since we can only measure the halo outside 
the point-spread function (at > 0f5 for our XMM-Newton 
data) , well beyond where dust associated with the source 
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would contribute. Below, we will find that the i ntensities 
of the halos confirm the extinction estimates of lDvK06bl 
and that, therefore, there is little extinction local to the 
source. In principle, it should also be possible to use 
the detailed structure of the halo to check the run of ex- 
tinction with distance inferred from the red-clump stars. 
Here, however, we find that, generally, the X-ray halos 
are fit just as well, or even better, with dust distribu- 
tions more uniform than indicated by the data, perhaps 
because simplifications made in our modelling of the dust 
removed some of the variety that would have smoothed 
the structure further. 

The layout of this paper is as follows. In §[51 we first de- 
scribe the model we use for fitting the X-ray observations. 
Next, in §[3j we describe how we measured the radial sur- 
face brightness profiles from archival XMM-Newton data 
and, in §[4j how we fit these with our model, discussing 
some consistency checks and possible pitfalls. We con- 
clude in §[5] 

2. EXPECTED SURFACE BRIGHTNESS DISTRIBUTIONS 

We model the observed surface brightness distributions 
for each source with three components, an unscattered 
central point source, a halo of light scattered by dust, 
and unrelated background emission. Below, we discuss 
our assumptions for the halo in more detail; for the 
point source, we will use measurements of an unabsorbed 
source, and for the background, we will assume it is in- 
dependent of position. 

The scattering of X rays on dust particles is, gener- 
ally, a complicated problem, and below we follow previ- 
ous work in making a number of simplifications, most of 
which are reasonable as long as one considers X rays with 
energies above ~ 1 keV. For our sources, this also implies 
the single-scattering limit is appropriate. 

Assuming that the dust grains are spherical, one can 
write the scattered intensity observed at an angular dis- 
tance 9 a hs from the central source as ()Smith fc Dwekl 
1998). 



Jobs 



F x / S{E)dE 



Nnn(a)da 



1 f(x) riff(F, q,6> scat ) 



dx, 



(1) 



where Fx and S(E) are the observed flux and normalised 
spectrum in the energy range [E m - m , F max ]; a is the grain 
radius and n(a) the corresponding number density; x is 
the distance from the observer normalis ed by the dis- 
tance to the source (see, e.g., Fig. 1 in iMathis fc Led 

1991); f(x) is the normalized spatial distribution of the 
scattering sites (equal to 1 if evenly distributed); and 
da (E , a, 9 sca t) / d£l is the differential scattering cross sec- 
tion, with S cat the angle over which the X ray is scat- 
tered. For small angles, 9 scat = #obs/(l — %)', hence, 
generally, the innermost regions of the halo are due to 
dust grains close to the source, while the outer regions 
are due to grains closer to the observer. 

In the above integral, the least-known parts are the 
dust size distribution and the scattering cross section, 
both of which depend on the physical state and com- 
position of the dust grains. For t he former, following 
IMathis. Rumpl. fc~N ordsicckl (|1977| ). generally a power- 



law form is assumed, while for the latter either the Gaus- 
sian approximation of the Rayleigh -Gans theory is used 
(hereafter RG; e.g.. I0verbecklll965f ) . or a more detaile d 
calculation is done using Mie theory (jvan de Hulstl ll957'). 
In the RG approximation, it is assumed that reflection 
from the grain surface is negligible and that the phase 
of the incident wave is not shifted inside the medium, 
resulting in co h erent addition of the scattered waves. 
iSmith fc Dwekl (jl998l ) have shown that for F > 1 keV 
and small scattering angles, the approximation agrees 
well with Mie theory, while for lower energies, RG the- 
ory underpredicts the scattering efficiency. 

In our work, we will assume the RG approximation, 
and limit ourselves to relatively high ener gies. Specifi- 
cally, we approximate the cross section by (|van de Hulstl 

Ell, 



da(E, o, 6> scat ) 

dn 



Cdusta exp 



scat 



(2) 



where the constant of proportionality Cdust depends on 
dust parameters such as the atomic c harge, mass number , 
density, and scattering factor (e.g... iHenke et al.l 11982ft . 
and <7 indicates the typical size of the halo produced 
by gr ains with radius a, given by ()Mauche fc Gorensteinl 

mm, 



a = 10.4 arcmin 



1 



(3) 



(F/lkeV)(a/0.lMm)' 

As an aside, we note that with this cross-section, and 
for f(x) — 1, the integral over x in Eq.[T]can be expressed 
in terms of the error function. Furthermore, the integral 
over annuli required to obtain mean surface brightnesses 
around some 9 oos , can also be expressed easily. Specifi- 
cally, 
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exp 
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dx 



cr 2 (l — x) exp(— a 2 ), 



(4) 



where we defined a = 9/^/2a(l — x). With this ex- 
pression, to calculate the expected halo brightnesses for 
piece- wise uniform dust distributions fix) , we only need 
to integrate numerically over a and F0 

For the dust size distribution, we use a power law 
n(a) oc a q with two different choices for the power and 
lower and upper limits. For our p rimary model , we se t 
q = -3.5, following the finding of IMathis et all (|1977fl . 
and use a mm = 0.00 5 um and tt ma = 0.25 / xm (informed 
by the work of, e.g.. iCostantini et al.ll2005l ). This model 
yields a halo domi nated by somewhat larger grains than 
found in the fits of lPredehl fc Schmittl (|1995h , whose cal- 
ibration of scattering optical depth with optical extinc- 
tion we will use. Therefore, to see how strongly our 
results depend on our choice, we also try a somewhat 
steeper power, q — —4, and smaller maximum grain size, 
flmax = 0.18 /mi, more simi lar to the values typically 
found by iPredehl fc SchmittJ . Finally, we fit profiles in 

1 For power-law dust distributions n(a) oc a q with integral power 
q, the integral over a can also be expressed in terms of error func- 
tions. 



3 



TABLE 1 
XMM-Newton Imaging Data Used 









Exp. 


b u 


Ay 


d 


Source 


Obs. ID 


Filter 


(ks) 


(°) 


(mag) 


(kpc) 


Anomalous X-ray Pulsars 














4U 0142+61 


0112781101 


Thin 


4.2 


-0.43 


3.5 + 0.4 


3.6 + 0.4 


IE 1048.1-5937 


0510010601 


Thin 


32 


-0.52 


5.6 : 


9.0 + 1.7 


1RXS J170849.0-400910 


0148690101 


Medium 


30 


+0.03 


7.7 + 2.2 


3.8 + 0.5 


Empirical PSF 














PKS 2155-304 


0080940301 


Thin 


39 


-52.25 








0158961301 


Medium 


39 


-52.25 







Note. — All observations used EPIC-pn in the small window mode. The exposure time is 
the total time not affected by flaring. For refer ence, the galactic latitude 6 n is l isted, as well 
as the extinction Ay and distance d inferred by Durant & van Kcrkwiik (2006a b). No error 
is listed for the extinction for IE 1048.1—5937, since it is based on a fit to broad-band X-ray 
data and the uncertainty is thus dominated by the extent to which the assumed intrinsic 
emission model is correct. 



narrow energy ranges, so that we can ignore the depen- 
dence on the spectral shape S(E) in our integration over 
energy. 

For given f(x), the above suffices to describe the radial 
profile of the halo, and our fits yield the normalisation. 
Integrating the shape over 8 b s then yields the total ob- 
served intensity of the halo, /halo which is related to the 
scattering optical depth r scat by, 



^halo 



/halo 



1 - exp(-T scat ), 



(5) 



where / core is the unscattered intensity of the point 
source. 

We normalize the scattering optical depth to 1 keV us- 
ing that r scat oc E~ 2 , and convert to equivalent visual 
extinction using the e mpirical calibration"! deriv ed from 
29 Galactic sources bv lPredehl fc Schmitti (fl995l) . 

7"scat , 1 kc V 

= A v (0.079 ± 0.003) - (0.052 ± 0.019). (6) 

We note that no offset would be expected, and its pres- 
ence may reflect the use of the RG approximation even 
at energies below 1 keV, where it is known to be inac- 
curate (P. Predehl, 2008, pers. comm.). The resulting 
systematic uncertainties are compensated to some ex- 
tent, however, by the fact th at the assumptions we mak e 
are very similar to those of iPredehl fc Schmittl (|1995l) . 
causing us to make errors in the same direction. Fur- 
thermore, assuming that the dust properties do not vary 
wildly between different lines of sight, the same should 
hold between our different sources, i.e., the relative in- 
tensities of the halos should provide a good measure of 
the relative dust columns. 

3. OBSERVED SURFACE BRIGHTNESS DISTRIBUTIONS 

In order to measure the scattering halos, we need good 
angular resolution and sensitivity, yet avoid non-linearity 
in the bright, unscattered light. These constraints are 
best met by data taken with the European Photon Imag- 
ing C amera with PN detectors (EPIC-PN; iStriider et all 
120011) onboard XMM-Newton (XMM). Most AXPs - and 
most good reference sources - are bright and thus ob- 
served with small-window mode to avoid pile-up. This 

2 The relation given here supercedes the different n umbers given 
in the text and abstract of lPredehl fc Schm itt (1995); P. Predehl, 
2008, pers. comm. 



helps to define the surface brightness accurately close 
to the source, but also implies we cannot measure the 
halo's full extent. In the XMM archive, suitable observa- 
tions are available for five of the six well-studied anoma- 
lous X-ray pulsars^ but two of those (IE 1841-045 and 
IE 2259+586) are associated with supernova remnants, 
the emission of which makes it very difficult to measure 
the scattering halo. Hence, we limited ourselves to the 
three remaining sources, IE 1048.1-5937, 4U 0142+61 
and 1RXS J170849.0-400910; see Table [Q As a refer- 
ence source to define the point-spread function, we used 
the BL Lac PKS 2155—304, which is at high galactic lat- 
itude and thus has a negligibly small scattering halo, and 
which has been observed regularly in different filters. 

For all sources, we reprocessed the raw data with the 
standard pipelines epchain and emchain from xmmsas 
version 7.1.0 (and the calibration files available with that 
release). We obtained the surface brightness distribu- 
tions by summing counts in concentric annular regions of 
10" width around the source, with the innermost annulus 
starting at 6'.'62 (avoiding the central parts, were system- 
atic differences due to small centering errors or pile-up 
effects (see below) are expected to dominate) and the 
outermost ending at 196'.'62. We selected events within 
given energy ranges using standard constraints for EPIC- 
PN (singles and doubles, and no warning flags). We cor- 
rected for out-of-time (OoT) events by subtracting an 
appropriately scaled radial profile, produced from simu- 
lated out-of-time only event lists (made with epchain; we 
used the scale factor of 0.011 listed in the XMM-Newton 
Users Handbook for small-window mode). We converted 
to surface brightnesses by dividing by the areas, using the 
arf gen and backscale tasks to ensure the areas include 
only the parts of the annuli that overlap the detector. 

A possible issue with our method, pointed out by the 
referee, is that the correction for out-of-time events might 
be unreliable if the source suffers from pile-up. This 
is because our analysis relies on the fact that photons 
detected "out of time" - during the transfer from the 

3 The sixth, XTE J1810-197, has only large-window data, for 
which there are no corresponding data of our point-spread function 
reference. We tried to analyze it using the PSF derived from small- 
window data, but could not find a good model for the scattering 
halo. This might be because of the different modes, the source's 
variability, or a break-down in our model. We did not investigate 
this further. 
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imaging part of the CCD to the storage area - have the 
same rate and spectrum as those measured during the 
integration. If pile-up occurs, however, the measured 
source event rate is reduced and the measured spectrum 
hardened, and thus one would underestimate the rate 
of low-energy, out-of-time events. Fortunately, for our 
observations, all taken in small-window mode, pile-up is 
not an issue (for our brightest source, 4U 0142+61, the 
total count rate is ^SOs" 1 , while in small window mode, 
with its read-out rate of 175 frames per second, pile- up 
becomes significant only for rates in excess of ~100s _1 ; 
we confirmed pile-up is not an issue using the epatplot 
task). Nevertheless, it would be a problem for brighter 
sources and/or imaging modes with longer integration 
times. In those cases, it would likely be best to measure 
the radial profile in areas excluding all parts affected by 
out-of-time events (i.e., excluding a strip along the read- 
out direction that passes over the source). 

4. SCATTERING HALOS 

To derive halo intensities, we fit our radial profiles with 
a combination of an empirical PSF, the scattering halo 
model for a variety of assumed distributions of dust along 
the line of sight (and for two dust size distributions; 
see §[2]), and a constant background. We fit the radial 
profile in two energy ranges simultaneously, 1.2-1.6 and 
2.0-2.4 keV. Here, the lower range will have the stronger 
dust scattering signal (yet is at sufficiently high energy 
that multiple scattering is not a significant issue for the 
sources considered here), but, given the small field of 
view, it is mostly sensitive to dust close to the source. 
The higher energy range will have a weaker scattering 
halo (though still noticeable), but helps to constrain the 
dust somewhat closer to the observer, as well as to verify 
consistency. 

We used the radial profiles from the observations of the 
high latitude source PKS 2155—304 as the empirical PSF 
in the fitting routines. The main advantages over using 
a theoretical PSF are that one automatically minimizes 
any common artefacts (cither instrumental or due to the 
analysis method; see, e.g., the wiggles at large radial 
distances in the right-hand panels of Figs [3J HI and 0, 
artefacts not necessarily captured in a theoretical PSF 
model. A disadvantage is that, effectively, the fits now 
measure relative brightnesses and background levels. In 
order to put these on an absolute scale, we determined 
the backgrounds and normalizations for our ca librator by 
fitting its profiles using the theoretical PSF of IGhizzardil 
(|2002t a King profile with energy and position-dependent 
width and exponent). 

For the dust distributions along the line of sight, we try 
two types. First, we use a uniform distribution, in which 
the halo contribution depends only on the total amount 
of dust. Second, we try a variety of piece- wise uniform 
distributions, constrained to follow the run of reddening 
with distance inferred from red clump stars bv lDv K06a. 
Here, the halo intensity depends on distance, which de- 
termines both the total dust column, and the fractional 
distances of the boundaries between the different pieces. 
Since the halo is fairly sensitive to the precise spatial 
distribution of dust, we chose up to three different dis- 
tributions. Fbr_the first - model A - we always tried to 
match the DvK06a results exactly (to the extent possible 
without having unphysical decreases of dust column with 
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Fig. 1. — Extinction and distance for 4U 0142+61. In the main 
panel, the run of extinc tion with distance is shown, as inferred 
from red-clump stars by |DvK06a l. Overdrawn are two piece-wise 
uniform approximations we used to fit the halo intensity (models 
A and B). In the top panels, we show for both models how the 
reduced x 2 of the fit to the observed radial profiles varies with 
distance, for two choices of dust distribution (as indicated). In all 
cases, two minima appear, between 4 and 6 kpc. In the right-hand 
panels, we show the implied variations in X 2 e( j with total extinction 
for both models and both dust distributions, as well as the results 
for a uniform dust distribution (for which the predicted scattering 
halos does not depend on distance). For all models, one infers 
A v ~ 4 mag. 




Fig. 2.— Radial profiles of 4U 0142+61 in the 1.2-1.6 keV (left) 
and 2.0-2.4 keV (right) ranges, fitted with a combination of an 
empirical PSF, a scattering halo for a variety of models, and a 
constant background. The green lines reflect the fit for a uniform 
distribution, with the dotted lines showing the unscattered and 
background components, and the solid lines the contribution from 
the halo and the total of all contributions. The best-fit halo con- 
tributions from other models are also shown (as labelled). In the 
two lower panels, deviations from the various fits are shown, with 
the residuals for fits using no halo or a halo for uniform dust in the 
middle panel, and those for fits for piece- wise linear dust distribu- 
tions in the lowest panel. Note that model A2 refers to the second 
minimum at large distance found for model A (see Fig. [Tt. 



distance). Since we found that sharp changes provided 
worse fits, we also tried up to two further models - B and 
C - both chosen to be "softer ," yet still consistent with 
the measurements of lDvK06al Combined, these models 
give a good sense of the sensitivity of our results to such 
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choices. 

In total, any individual fit has five free parameters: 
two unscattered source intensities and two background 
brightnesses (one for each of the two energy ranges), and 
either the total dust column or distance (which deter- 
mines the fractional halo intensities). In practice, in our 
fits for uniform density, the five parameters are deter- 
mined simultaneously by linear least-squares, while for 
those with varying dust distributions, we step through 
distance, and for each distance calculate the expected 
fractional halo contribution and fit for the corresponding 
best-fit unscattered source intensities and backgrounds. 

We show our fit results in Figs. [THB1 and summarize 
numerical results for our primary dust distribution (with 
q = —3.5 and a max = 0.25 fim) in Table [2j Two conclu- 
sions stand out: (i) fits without a halo are much poorer 
than those with a halo; and (ii) the fractional intensity 
of the halo increases as expected from the reddening in- 
ferred from the X-ray spectra, with 4U 0142+61 having 
the smallest halo signal, IE 1048.1—5937 an intermediate 
one, and 1RXS J170849.0-400910 the largest. Below, we 
describe the results for the individual sources in detail, 
and then discuss possible model caveats and systematic 
uncertainties. 

4.1. 4U 0142+61 

One of the brightest and most studied AXPs is 4U 
0142+61. Absorption edges in its high-resolution X- 
ray spectr um indicate a total extinction Ay = 3.5 ± 0.4 
(|DvK06bf ). which, combined with the run of extinction 
with distance inferred from the red clump stars, suggests 
a distance of about 3.5 kpc, in a region of rapidl y increas- 
ing optical extinction, presumably a spiral arm (|DvK06at 
see Fig. |T]). Thus, we expect that a large fraction of the 
dust along the line of sight is close to the source, which 
should affect the inner parts of the halo in particular. 

In Fig. [TJ we show how \ 2 varies as a function of 
reddening and distance. For uniformly distributed dust 
with our primary dust size distribution (q = —3.5 and 
imax = 0.25 fJxa), there is a clear minimum around 
Ay — 4.2 mag, close to the value inferred spectroscop- 
ically. The piecewise uniform models also show consis- 
tent minima, with reddening of 3-4 mag and distances 
of ^3.6 kpc. However, they also have secondary minima 
at larger distances, of ~ 6.7 kpc. Given that AXPs are 
young and thus likely to reside in spiral arms, the first 
minimum almost certainly is the correct one. 

For our second dust size distribution (q = —4 and 
Umax = 0.18 /im), we find fits of similar quality, but some- 
what larger values of the extinction. This is as expected: 
for smaller grains, more of the halo will be outside our 
field of view, and hence a stronger halo is required to 
match the observed halo contribution inside the field of 
view. From these, we conclude the uncertainty due to the 
dust distribution is considerable, although the results do 
not change qualitatively. 

In Fig. [2 we show the fits to the radial profile, as well 
as the residuals. From the different halo curves, one sees 
how larger quantities of dust concentrated towards the 
source increase the innermost parts of the halo compared 
to the outer wings. Indeed, both piecewise models give 
very poor fits if the source is placed near the far end or 
just behind the rapid rise of extinction around 3.7 kpc 
inferred from the red clump stars. For the best fits, how- 
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Fig. 3.— Like Fig. [T] but for IE 1048.1-5937. Note that both 
approximations to the run of reddening with distance use linear 
extrapolation beyond Ay ~ 6 mag. The fits indicate a total ex- 
tinction of about 6 mag for our primary dust size distribution, and 
a distance of ~ 10 kpc. 



ever, the residuals are rather similar, since the changes 
in halo shape can be compensated to some extent by 
changes in the other parameters, in particular the back- 
ground. Note that in the higher energy range, the halo is 
already very weak compared to the source: its fractional 
intensity is about 5%. 

One surprise is that the uniform model gives a better 
fit than the two piecewise uniform ones, even though 
the latter presumably are closer to reality (as inferred 
from the red-clump stars). Below, we find that smoother 
distributions also fit better for the other sources. We 
briefly return to this in §[5j but stress here again that 
among all models, the total extinction is fairly similar. 

4.2. IE 1048.1-5937 

One of the main goals of the present study was 
to test the extinction and distance estimates for IE 
1048.1—5937. Based on the coincidence with a H I bub- 
ble, Gaensler et al. suggested a distance of ~ 2.7 kpc, but 
the high extinction measured for this source, Ay ~ 6, 
combined with the run of extinction wi th distanc e, sug- 
gests a much larger distance of ~ 9 kpc (|DvK06a| ) . 

A major uncertainty is the estimate of the extinction, 
since it was based on broad-band X-ray fits, and thus 
necessarily dependent on the assumed underlying source 
model. Unfortunately, the source was not bright enough 
to detect absorption edges in its grating X-ray spectra 
(DvKOGb). Visual inspection of observed radial profile 
of IE 1048.1—5937, however, already suggests the ex- 
tinction is indeed fairly high: the source clearly has a 
contribution due to a scattering halo at least as large as 
4U 0142+61 (compare, e.g., the offset around 60" from 
the point-spread function in Figs [2] and . 

This is confirmed by our scattering halo fits (Fig. [3|): 
all require Ay around 6 mag (somewhat higher for our 
secondary choice of dust distribution; note, though, that 
for those we had to extrapolate the run of extinction 
with distance beyond the range covered by the red clump 
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Fig. 4.— Like Fig. [2] but for IE 1048.1-5937. 

stars). Thus, like for 4U 0142+61, we find that the scat- 
tering halos are consistent with the total extinction in- 
ferred from the X-ray spectrum, and, based on the red- 
clump stars, the inferred distance remains high, ~ 10 kpc. 

The high extinction (and strong halo) are inconsistent 
with the extinction found from the red-clump stars at 
2.7 kpc. The explanation can also not be that there 
is strong extinction close to the source, since this dust 
would contribute only to the very inner parts of the halo, 
and thus not to the halo measured in our fit. 

Comparing the different halo models in detail, we find 
again that the fit for the more uniform dust distributions 
is the best. Indeed, following precisely t he run of extinc- 
tion with distance inferred bv lDvK06al . the fit is rather 
poor, with reduced x 2 — 1-7. From the halos and residu- 
als in Fig. 21 one again sees that the sharper increase in 
Ay near the source leads to a halo that is less consistent 
with the observed one at small angles. 

4.3. 1RXS J170849.0-400910 

Along the line of sight to 1RXS J170849.0-400910, 
IDvK06al inferred large jumps in reddening from the red 
clump stars (Fig. [5]), likely associated with the Carina 
and Crux spiral arms. Given Ay = 7.7 ± 2 .2 inferre d 
from absorption edges in the X-ray spectrum (DvK06b), 
the source should be inside the Crux spiral arm, at 3.8 ± 
0.5 kpc. 

Our fits generally reproduce the above, with Ay be- 
tween 8 and 9 mag, and the distance around 4 kpc. 
Again, however, the fits are much better for smoother 
dust distributions (Fig. [SJ Table [2]). Indeed, if we as- 
sume the exact dust distribution inferred from the red 
clump stars, we obtain an unacceptable fit, with reduced 
X 2 — 2.6. In Fig. [SJ one sees that the poor fit is due to 
the sharp increase in reddening near the source, which 
causes a halo that is too strong and too flat at small an- 
gles compared to what is required by the observations. 

4.4. Background Consistency Check 

Within the small, ~ 256 x 260 arcsec area covered by 
our small-window mode observations, the profiles barely 
reach the background, even without a scattering halo. 
With a halo, which has a typical size 2a ~ 1000" at 
1.4 keV, the profile extends well outside. Thus, one can- 
not determine the background from an "empty" region, 
which is why we left it as a free parameter. 

Nevertheless, some consistency checks are possible. 
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Easiest perhaps would be to compare with background 
estimates found in the outer regions of EPIC-MOS 
exposures taken simultaneously with the PN observa- 
tions. Unfortunately, we found that those ratios were 
not constant, probably at least in part due to detec- 
tor and position dep endent instrumental backgrounds 
(jCarter fc Read! 120071 ) . A rough sense can be obtained 
from blank fields, even though no templates are avail- 
able for small-window mode. These show that the typi- 
cal backgrounds in the two energy bands considered are 
both ~ 0.5 x 10~ 6 s~ x arcsec -2 , with at least a factor 
2 variation between different observations These are 
roughly consistent with the rates we infer from our fits 
(see Table EJ. 

A stronger constraint might be set by comparing the 
rates betwe en different ener g y ban ds. This is because, 
as found bv lKatavama et aT] (|2004fl . the spectrum of the 
background has no strong dependence on intensity, at 
least at > 1.7 keV. Indeed, for our calibrator, we find that 



* http : //xmm2 . esac . esa. int/ext ernal/xmm_sw_cal/background/ index. shtnL 
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X-ray Scattering Halo Fits 
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Note. — These fits are for a dust size distribution n(a) <x a~ 3,5 
between a m i n = 0.005 and a max = 0.25 fim. We omitted formal uncer- 
tainties, since they are typically substantially smaller than the differ- 
ences between different models. The fluxes / are those of the unscat- 
tered component; the background surface brightnesses b have units of 
10 — 6 s — 1 arcsec -2 . 

a For the 4U 0142+61 models with a '2' suffix, the solution is that of 
a second minimum. See Fig. [2] 

b The number of degrees of freedom is 34 for the fits without halo, and 
33 for fits with a halo model (38 data points, 4 or 5 fit parameters). 



we reproduce the ratio of ~ 1.1 found bv lKatavama et all 
between the background rate in the 4-7.2 keV range and 
that in the 1.7-4 keV range. From our calibrator, we find 
that for the bands we use, the ratio between the back- 
ground counts with those in the 4-7.2 keV band are ^0.4 
and 0.2 (for 1.2-1.6 and 2.0-2.4keV, respectively; the 
ratio for the lower energy band is more un c ertain , con- 
sistent with the finding of iKatavama et aT] ([2004), who 
found that the background rate below 1.7 keV was poorly 
correlated with the rate at higher energies) . We also find 
that all observations have similar 4-7.2 keV background 
count rates, of 1.5-2.5 x 10 -6 s -1 arcsec -2 . Thus, the ex- 
pected background rates in the 1.2-1.6 and 2.0-2.4 bands 
are also similar, ~0.8 and 0.4 x 10 -6 s _1 arcsec -2 . This 
is roughly consistent both to what we listed above for 
the blank fields, and with what we infer from those of 
our fits that had acceptable \ 2 - 

Overall, we conclude that the fact that we cannot mea- 
sure the background rate directly, and thus have to leave 
it as a free parameter in our fits, affects the precision of 
our analysis, but has not introduced any major system- 
atic uncertainties. 

4.5. Model Caveats & Uncertainties 

Our approach to fitting the halos has a number of po- 
tential pitfalls. First, our observations cover only part 
of the scattering halo: within the ~ 3' covered by the 
data, only about 60% of the halo signal is contained (for 
uniformly distributed dust). However, to first order, this 
should just lead to large uncertainties, not large system- 
atic errors: since we leave the background free, our lim- 
ited coverage simply means it is easier to obtain an ade- 
quate fit. Furthermore, for our higher energy band, a bit 



more of the halo is captured (about 75% for uniformly 
distributed dust). 

Second, the model we used to fit the observed scatter- 
ing halo is rather simplified and thus it is worthwhile to 
see what the possible systematic errors could be. One 
relatively small simplification is that we take the mean 
energy of all our photons to be that of the middle of a 
given energy band (and also use this to scale to the scat- 
tering optical depth from that at 1 keV) . The associated 
error is likely small, however, since the count spectra of 
our sources show that the photon energy distribution is 
quite uniform over our 0.4 keV bands. (Note that we 
do integrate over the band for calculating the scattering 
halo, so include the small amount of "smearing" because 
of the change of energy.) 

Another approxim ation we make is th at we ignore 
multiple scattering. iMathis fc Led (|1991f ) have shown 
that this dominates for r scat > 1.3 or Ay > 17 mag 
at IkeV. For our lowest energy band, at 1.4 keV, the 
use of the single-scattering approximation thus requires 
Ay <C 33 mag, which is reasonably well satisfied for our 
sources. 

The largest simplifications we make are the assumption 
of spherical dust grains with a power-law size distribu- 
tion, and the treatment of scattering in the RG approxi- 
mation. For the dust distribution, our two trial assump- 
tions give significantly different results. We tried leav- 
ing the power-law index and limits free in our fits, but 
did not reach a clear conclusion. Generally, though, we 
would expect any systematic errors should be in the same 
direction for all sources, (unless dust properties vary sig- 
nificantly for different lines of sight, e.g., because one 
crosses a molecular cloud while another does not). Thus, 



e.g., our conclusion that the total dust column towards 
IE 1048.1-5937 is about 50% larger than that towards 
4U 0142+61 should be reasonably safe. Indeed, as al- 
ready mentioned in §[2j since we us e the calibration of 
the ha lo intensity with extinction of IPredehl fc Schmittl 
(1995), systematic errors may well be minimized by us- 
ing simplifications as similar to theirs as possible, as we 
do. 

Despite the above, ideally one would use phys- 
ical inodel_s_^oj__yie__dust ! __such as that proposed 
by iWeingartner fc Draind ()2001l ). Unfortunately, at 
present, our current understanding o f dust scatter ing ap- 
pears to be inadequate. For instance. ISmithj ( [20081 ) found 
that no model could reproduce in detail the dust scatter- 
ing halo around GX 13+1, which he c arefully measured 
in annuli ranging from 2 to 1000 arcsec. iSmithl also found 
that uniform dust distributions with power-law grain size 
distributions fitted the data as well as the presumably 
more physical models. Thus, at present, it seems one 
will have to continue to rely on empirical calibrations 
between halo intensity and optical depth or reddening. 

5. CONCLUSIONS 

We used XMM observations to measure the strengths 
of the scattering halos for three Anomalous X-ray Pul- 
sars. In all three cases, we find that they are consistent 
with what was expected given extinctions determined 
previously from fits to the X-ray spectra. Combined with 
the run of extinction wi th distan ce, we thus also confirm 
the distances found by DvK06a. In particular, our re- 
sults confirm that IE 1048.1—5937 is much further than 



inferred bv lGaensler et a l. (2005) from a possible associ- 
ation with an H I bubble. 

Somewhat surprising is that the observed radial pro- 
files are consistently fit as well or better by uniform dust 
distributions than they arc by dust distributions that try 
to follow the actual run of extinction with distance in- 
ferred from red-clump stars. Part of the reason for the 
poorer fits for the presumably more realistic dust distri- 
butions may be that in our analysis we ignored complex- 
ities in the dust properties or the scattering process that 
might tend to smooth the radial profile, mimicking what 
would result from a more uniform dust distribution along 
the line of sight. 

Nevertheless, our results suggest that, if one calibrates 
one's analysis method with a source with known extinc- 
tion (4U 0142+61 in our case), it is possible to infer 
estimates of extinction columns to other sources good 
to about 20% from their scattering halos. With proper 
calibration, it should thus be possible to use the halos 
as independent constraints on model fits to the X-ray 
spectrum, and to help determine distances by comparing 
with the run of extinction with distance. 
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